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ABSTRACT 

We present the measurements and modelling of the small-to-intermediate scale 0.1- 
25 (i“^Mpc) projected and three-dimensional (3D) redshift-space two-point correlation 
functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data 
Release 7. We find a clear dependence of galaxy clustering on luminosity in both 
projected and redshift spaces, generally being stronger for more luminous samples. 
The measurements are successfully interpreted within the halo occupation distribu¬ 
tion (HOD) framework with central and satellite velocity bias parameters to describe 
galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. 
In agreement with previous studies, we find that more luminous galaxies reside in 
more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD con¬ 
straints. Moreover, we find that luminous central galaxies are not at rest at the halo 
centres, with the velocity dispersion about 30% that of the dark matter. Such a rela¬ 
tive motion may reflect the consequence of galaxy and halo mergers, and we find that 
central galaxies in lower mass haloes tend to be more relaxed with respect to their host 
haloes. The motion of satellite galaxies in luminous samples is consistent with their 
following that of the dark matter. For faint samples, satellites tends to have slower 
motion, with velocity dispersion inside haloes about 85% that of the dark matter. 
We discuss possible applications of the velocity bias constraints on studying galaxy 
evolution and cosmology. In the appendix, we characterize the distribution of galaxy 
redshift measurement errors, which is well described by a Gaussian-convolved double 
exponential distribution. 
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1 INTRODUCTION 

The three-dimensional galaxy distribution in our Universe 
can be probed through the large-scale galaxy redshift 
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surveys, such as the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). The angular positions of the galaxies can 
be accurately measured in photometric observations, while 
the radial positions are usually obtained from the observed 
galaxy redshifts. However, the radial distances derived from 
the redshifts differ from the real positions of galaxies due 
to the existence of the galaxy peculiar velocities, which is 
usually referred to as the redshift-space distortion (RSD) 
effect. Although the RSD prevents us from measuring the 
true galaxy distribution, it also provides valuable informa¬ 
tion about the galaxy kinematics in dark matter haloes. 

A commonly used statistic for analysing the galaxy dis¬ 
tribution is the three-dimensional (3D) two-point correlation 
function (2PCF), ^(r'p,r.,r), where Vp and r-^ are the trans¬ 
verse and line-of-sight (LOS) separations of galaxy pairs, re¬ 
spectively (see e.g. Zehavi et al. 2005, 2011; Li et al. 2006; 
Guo et al. 2013). To minimise the effect of RSD, the tradi¬ 
tional way of probing the real-space galaxy distribution is 
through the projected 2PCF, Wp{rp), by integrating ^{vp, r^) 
along the LOS direction. Since Wpirp) is insensitive to the 
galaxy peculiar velocities, it only probes the galaxy spatial 
distribution. To fully understand the galaxy phase-space dis¬ 
tribution, we need to model the redshift-space clustering of 
galaxies. 

In contemporary galaxy formation and evolution mod¬ 
els, galaxies form and evolve in dark matter haloes. The 
galaxy distribution in the universe can then be studied with 
the halo model, through the distribution of galaxies within 
the haloes and the halo distributions in the large-scale struc¬ 
ture of the universe (see e.g. Cooray & Sheth 2002, and 
references therein). Since the distribution of dark matter 
haloes is well understood using analytic models and nu¬ 
merical A^-body simulations (e.g.. Mo et al. 1996; Springel 
2005; Klypin et al. 2014), the key component of the mod¬ 
els is the connection between galaxies and dark matter 
haloes, such as the framework of halo occupation distri¬ 
bution (HOD) or the closely related conditional luminos¬ 
ity function (e.g. Jing et al. 1998; Peacock & Smith 2000; 
Seljak 2000; Scoccimarro et al. 2001; Berlind & Weinberg 
2002; Yang et al. 2003; Zheng et al. 2005; Guo et al. 2014; 
Skibba et al. 2015). The HOD describes the probability dis¬ 
tribution P{N\M) of having N galaxies of a given type in 
a dark matter halo of virial mass M. The probability dis¬ 
tribution P{N\M), together with the spatial and velocity 
distributions of galaxies inside haloes, is crucial to interpret 
and understand the real- and redshift-space clustering of 
galaxies. The observationally inferred HOD can help us test 
and constrain galaxy formation models. 

A large fraction of the information contents on the HOD 
(like the occupation function and galaxy kinematics and spa¬ 
tial distribution inside haloes) are contained in the small- 
scale clustering of galaxies. However, both measuring and 
modelling the galaxy distribution are nontrivial in the small- 
scale nonlinear regime. On the observational side, in fibre- 
fed spectrograph surveys as in the SDSS, the hardware limit 
that two fibres on the same plate cannot be placed closer 
than an angular separation of 55” (Blanton et al. 2003a) sig¬ 
nificantly hinders the number of close galaxy pairs on small 
scales. Fortunately, this fibre collision effect can be accu¬ 
rately corrected using the method of Guo, Zehavi & Zheng 
(2012), by taking advantage of the recovered redshifts of 
collided galaxies in the plate-overlap regions. On the theory 
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Figure 1. Construction of volume-limited galaxy samples with 
different luminosity thresholds. The dots show the distribution 
of SDSS galaxies as a function of redshift z and r-band absolute 
magnitude Mr- The different colour lines delineate the selection 
cuts for different luminosity-threshold samples defined in Table 1. 

side, modelling the galaxy distributions on small scales is 
also difficult, especially in redshift space due to the lack of 
understanding of the galaxy phase-space distribution. This 
problem can be alleviated with the help of high-resolution 
A-body dark matter simulations. 

Recently, Guo et al. (2015) (hereafter G15) measured 
and modelled the luminous red galaxy (LRG) distribution 
in the SDSS-HI Baryon Oscillation Spectroscopic Survey 
(BOSS; Eisenstein et al. 2011; Dawson et al. 2013) at red¬ 
shift z~0.5, and found that central galaxies are not at rest at 
the halo centres and the satellite galaxies move more slowly 
than the dark matter. The difference in galaxy and mat¬ 
ter velocity distributions is dubbed as velocity bias. In this 
paper, we follow the method of G15 and infer the phase- 
space distribution of galaxies in the local universe through 
modelling the redshift-space galaxy 2PCFs. We improve the 
model of G15 and also incorporate a more accurate red¬ 
shift error model shown in Appendix A. In Section 2, we de¬ 
scribe the data, the galaxy samples, and the redshift-space 
2PGF measurements. We introduce our modelling method 
in Section 3. The constraints on the occupation function and 
galaxy phase-space distribution are presented in Section 4. 
Finally, we summarize our results in Section 5. 

Throughout the paper, for the measurements we assume 
a spatially flat A cold dark matter (AGDM) cosmology, with 
Dm = 0.307 and h = 0.678, which is adopted in the Multi- 
Dark simulation (and model) we use, and consistent with the 
constraints from Planck (Planck Collaboration et al. 2014). 


2 DATA AND MEASUREMENTS 

For the purpose of studying the galaxy distribution in 
the local universe, we use the galaxy sample of the New 
York University Value-Added Galaxy Catalog (NYU-VAGC; 
Blanton et al. 2005), which is constructed from the SDSS 
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Data Release 7 Main galaxy sample (Abazajian et al. 2009). 
The sample covers an effective area of about 7300 deg^, with 
galaxies selected using an r-band Petrosian magnitude limit 
of r < 17.77. The magnitudes in the catalog are A-corrected 
and passively evolving to the median redshift of z = 0.1 
(Blanton et al. 2003b). To properly measure and model the 
galaxy clustering and its dependence on galaxy luminosity, 
we construct volume-limited samples with different luminos¬ 
ity thresholds. We impose a minimum redshift of z = 0.02. 
The selection cuts for the different samples are shown in 
Fig. 1. Table 1 provides the corresponding sample informa¬ 
tion, including the average number density and the volume 
of each sample. Overall, the samples are similar to those 
constructed in Zehavi et al. (2011). 

We measure the 3D redshift-space 2PCFs 5(rp,r,r) 
of the SDSS DR7 Main galaxy sample through the 
Landy-Szalay estimator (Tandy & Szalay 1993). The 2PCF, 
^(rp, Tt,), can be further integrated along the LOS to reduce 
the effect of RSD. The resulting projected 2PCF Wp{rp) 
(Davis & Peebles 1983) is defined as 

noo 

mp(rp) — ^ ^ . (1) 

Jo 

In both the measurements and the model, we integrate 
^(’■p!^ir) to r.,r,max = 40/i“^Mpc to obtain Wp. We adopt 
logarithmic Vp bins centred at 0.13 to 20.48/i“^Mpc with 
A log Tp = 0.2, and linear bins from 0 to 40 /i“^Mpc with 
Ar^ = 2 /i“^Mpc. Denoting the relative redshift-space posi¬ 
tion of a pair of galaxies as s, we also measure the redshift- 
space 2PCF in bins of s and jj., with s = being the 

redshift-space separation of galaxy pairs and fi being the co¬ 
sine of the angle between s and the LOS. The redshift-space 
2PCF ^(s, jj,) can be expanded into multipoles (Hamilton 
1992), 

( 2 ) 

i 

where Pi is the l-th order Legendre polynomial. The mul¬ 
tipole moments are usually used to characterize the 
redshift-space clustering (G15). We focus on the measure¬ 
ments of the monopole (^o), quadrupole (^ 2 ), and hexade- 
capole (^ 4 )- For s, we adopt the same logarithmic bins as rp, 
while for fj, we use linear bins from -1 to 1 with Ajj. = 0.05. 

In linear theory, the three multipole moments we adopt 
are the only non-zero terms. While at any scales odd terms 
are zero by the symmetry of s-t fhe translinear or 

nonlinear scales explored in this paper, higher-order even 
terms do exist. However, the information content in the 
higher-order terms is minimal compared to the three main 
terms (^o, C 2 , and ^ 4 ), as they are highly correlated with 
these lower-order terms. For example, Hikage (2014) ex¬ 
plored the constraints on the HOD by including multi¬ 
poles of different orders, and found that including the tetra- 
hexadecapole {I = 6 ) in addition to multipoles with I = 0 , 2 , 
and 4 leads to virtually no improvement in the constraints 
(see their Table 1 and Fig. 3). We therefore limit our study 
to only ^ 0 , C 2 , and ^ 4 , in addition to Wp. 

We apply the method developed in 
Guo, Zehavi & Zheng (2012) to correct the fibre-collision 
effect, enabling accurate measurements of the small-scale 
2PGFs. For each sample, the covariance matrix of the 


Table 1. Samples of different luminosity thresholds 


jy^max 

■^max 

N 

rig 

V 

-18.0 

0.041 

35649 

31.42 ±4.03 

1.19 

-18.5 

0.053 

57786 

22.25 ±2.70 

2.72 

-19.0 

0.064 

72484 

15.66 ±2.06 

4.87 

-19.5 

0.085 

125664 

11.64 ± 1.27 

11.44 

-20.0 

0.106 

131623 

6.37 ±0.75 

22.00 

-20.5 

0.132 

122115 

3.13 ±0.30 

41.82 

-21.0 

0.159 

76802 

1.16 ±0.12 

71.74 

-21.5 

0.198 

35207 

0.29 ± 0.03 

134.65 


The absolute magnitude is computed by assuming h = 1. The 
minimum redshift of all the luminosity threshold samples is 
■^min = 0.02. The total number of galaxies N in each sam¬ 
ple is also displayed. The mean number density rig is in units 
of 10“®/i®Mpc“®. The volume V of each sample is in units of 
lO^h-^Mpc®. 

measurements is estimated from 400 jackknife samples 
(Zehavi et al. 2011; Guo et al. 2015). 

We show in Fig. 2 the measured 2PGFs for the different 
luminosity threshold samples. It is clear from the projected 
2PCF Wp(rp) and the redshift-space monopole ^o(s) that 
more luminous galaxies have stronger clustering amplitudes 
than their fainter counterparts, consistent with the results of 
Zehavi et al. (2011). The measurements of the quadrupole 
^ 2 ( 5 ) and hexadecapole ^ 4 (s) are noisier, but the overall de¬ 
pendence on luminosity is clear. For example, more luminous 
galaxies show a more positive quadrupole and hexadecapole 
on small scales, indicating stronger Fingers-of-God (FOG) 
effects. We jointly model all the 2PCF measurements in the 
following sections. 


3 SIMULATION AND MODEL 

To accurately model the galaxy clustering under the HOD 
framework, we follow the simulation-based model method 
developed in Zheng & Guo (2015), which is used in G15 to 
model redshift-space clustering of BOSS galaxies. It tabu¬ 
lates properties of haloes in an A-body simulation (e.g., halo 
mass functions, halo profiles, halo clustering) necessary for 
computing galaxy 2PGFs. Given a set of HOD parameters, 
the model can then accurately predict the galaxy 2PGFs. It 
is equivalent to assigning galaxies to haloes in the simula¬ 
tion with the given set of HOD parameters and measuring 
the 2PGFs from the resultant mock catalog. When popu¬ 
lating haloes with galaxies, the redshift-space distortion is 
applied in a galaxy-by-galaxy manner by using the velocity 
information of galaxies and haloes. Our modelling method 
does the same thing by using the galaxy velocity distribu¬ 
tion inside haloes and the redshift-space clustering of haloes 
(see Zheng & Guo 2015 for details). Our simulation-based 
modelling method works more efficiently than directly pop¬ 
ulating the simulation with galaxies, since it enables ‘pop¬ 
ulating galaxies’ and ‘measuring the 2PGF in the mock’ to 
be performed analytically. The model is accurate since it au¬ 
tomatically takes into account the effects of halo exclusion, 
nonlinear growth, and scale-dependent halo bias by using 
the halo catalogues in high-resolution simulations. In par¬ 
ticular, it is well suited to model the redshift-space galaxy 
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Figure 2. Measured projected 2PCF Wp{rp) (top left), redshift-space monopole (top right), quadrupole (bottom left), and 

hexadecapole ^ 4 (s) (bottom right) for SDSS galaxies. The different colour lines show the measurements for different luminosity threshold 
samples. The errors are estimated from 400 jackknife samples. 


clustering on small and intermediate scales, for which an 
accurate analytic model is difficult to develop (e.g. Tinker 
2007; Reid & White 2011; Wang et al. 2014a; White et al. 
2015). 

The model we use in this paper is based on the Multi- 
Dark simulation of Planck cosmology (MDPL; Klypin et al. 
2014), which is carried out with L-GADGET-2 code 
(Springel 2005). The cosmological parameters (Dm = 0.307, 
D(, = 0.048, h = 0.678, Ua = 0.96, and as = 0.823) 
used in MDPL are consistent with the recent results from 
Planck (Planck Gollaboration et al. 2014). The simulation 
has 3840^ dark matter particles in a box of lh~^ Gpc 
(comoving) on a side, so the mass resolution is 1.51 x 
10 ® /i“^Mq, which is about 6 times higher than the previ¬ 
ous MultiDark run simulation (MDRl; Prada et al. 2012). 
The force resolution, i.e. the gravitational softening length, 
is only 5/i“^ kpc (physical) at low redshifts, which enables 
us to accurately model the clustering signals on very small 
scales. We use the simulation output at « = 0 to model all 
the luminosity threshold galaxy samples in the NYU-VAGC. 
In principle, when modelling the measurements, it is better 
to choose the simulation output to match the mean redshift 
for each individual sample. This is certainly limited by the 
available simulation outputs. On the other hand, given the 
small redshift range of the SDSS Main galaxies, the effect of 


using one output (as we do) is small. We have tested apply¬ 
ing the z = 0.1 simulation output for modelling the data, 
and the inferred HOD parameters are consistent with those 
from the default model built on the z = 0 output. 

Dark matter haloes are identified using the Rock- 
star phase-space halo finder (Behroozi et al. 2013), which 
is efficient and accurate in finding the bound spherical 
structures from the density peaks in the phase space 
(Onions et al. 2012; Knebe et al. 2013). Halo mass is defined 
from the given spherical overdensities of a virial structure 
(Bryan & Norman 1998). Note that we do not remove the 
unbound particles in the haloes, because the satellite galax¬ 
ies in the haloes can also be unbound. The halo positions, 
velocities and velocity dispersions are calculated from all the 
particles in the haloes. In Rockstar, the centre of each halo 
is computed from the average particle locations for the in¬ 
ner friends-of-friends (FOF) subgroup that best minimizes 
the Poisson error. Different from G15, who use the aver¬ 
age velocity of the inner 25 per cent of halo particles as 
the halo core velocity, we define the halo velocity as the 
average velocity of all particles in the halo, i.e. the centre- 
of-mass velocity. The purpose of this definition is to make 
better comparisons with the literature, and also make easier 
the application of our models to other low-resolution simu¬ 
lations. The definition of the halo velocity is important for 
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comparing the results of galaxy velocity bias, since the halo 
core velocity can have a substantial velocity offset from the 
halo bulk velocity (Behroozi et al. 2013; Reid et al. 2014). 

For the HOD modelling of the galaxy clustering, we 
follow the parametrization of Zheng, Coil & Zehavi (2007) 
by decomposiirg the contributions to the mean occupation 
function {N{M)) of galaxies (i.e. the average number N of 
galaxies in a sample in haloes of mass M) into the central 
and satellite components. 


{N{M)) 


(iVcen(M)) 


(Arcen(M)) + (iVsat(M)), (3) 

1 / log M - log Mmin 

2 _ \ 



(iVsat(M)} = (iV_(M)} (5) 

where Mmin describes the cutoff halo mass of the central 
galaxies and criogM takes into account the scatter between 
the galaxy luminosity and halo mass. The three parame¬ 
ters for the satellite galaxies are the cutoff mass scale Mo, 
the normalization mass scale M[ and the power-law slope 
a at the high-mass end. In our model, we implicitly assume 
that the halo hosting a satellite galaxy in a given luminos¬ 
ity threshold sample also hosts a central galaxy from the 
same sample. One derived parameter we have is Mi, the 
characteristic mass of haloes hosting on average one satel¬ 
lite galaxy. In combination with the halo mass function, the 
satellite fraction /sat of the galaxies in the sample can also 
be derived from the central and satellite mean occupation 
functions. 

To model the redshift-space galaxy clustering, we need 
to specify the phase-space (spatial and velocity) distribution 
of galaxies inside haloes. We put the central galaxies at the 
halo centres and randomly select dark matter particles inside 
haloes to represent the satellite galaxies. When calculating 
the redshift-space clustering, we employ the plane-parallel 
approximation and use the z direction in the simulation as 
the LOS. The shift of a galaxy’s position from real space 
to redshift space in the z direction due to the RSD is then 
calculated as AZ = 114 ( 1 - 1 - z)!H{z) with 2 = 0, where Vi is 
the LOS peculiar velocity of the galaxy. 

As found by G15, the motion of galaxies can differ from 
that of dark matter. We therefore introduce two velocity 
bias parameters (oc and Os) in the HOD model when de¬ 
scribing the velocities of central and satellite galaxies. Note 
that all the velocity quantities below refer to the 2 (LOS) 
component for our modelling purpose, including the central 
galaxy velocity Uc, satellite galaxy velocity Vs, halo velocity 
Vh, and particle velocity dispersion (t„ inside a given halo. 
We first measure the LOS velocity dispersion (t„ from all 
the dark matter particles in the haloes. The central galaxy 
is not necessarily at rest with respect to the host halo, and 
its velocity Vc — Vh in the frame of the halo is assumed to 
follow a Laplace distribution, in the form of 


f{vc - Vh) = exp 

V 



( 6 ) 


where Vh is the LOS centre-of-mass velocity of the halo, Uc = 
OcUv is the LOS central galaxy velocity dispersion, and Oc is 
the central galaxy velocity bias, characterizing the relative 
motion between the central galaxy and the host halo. The 
use of Laplace distribution instead of the commonly-used 
Gaussian distribution is motivated by the distribution of 


the velocities of brightest cluster galaxies (BCGs) relative 
to satellites in Abell clusters (Lauer et al. 2014). 

To allow for possible velocity offsets between the satel¬ 
lite galaxies and the randomly-selected dark matter parti¬ 
cles, we scale the velocity of the satellite galaxies in the 
centre-of-mass frame of the halo by a satellite velocity bias 
factor Os, 


Vs - Vh = OLs{Vp - Vh), (7) 

where Vs and Vp are the LOS velocities of the satellite 
galaxies and the selected dark matter particles, respectively. 
Therefore, the LOS velocity dispersion as of satellite galax¬ 
ies in the haloes is Os = Q-sOn (see, e.g. Tinker 2007). We 
note that even though we only apply the LOS velocity bias 
in the above equations, the velocity bias exists in all compo¬ 
nents of the galaxy velocities. However, for the purpose of 
modelling the redshift-space clustering, only the LOS com¬ 
ponent matters. In our fiducial HOD model, we assume a 
constant galaxy velocity bias, good enough given the cur¬ 
rent data precision. 

Except for the galaxy velocity bias, another ingredient 
that could affect the galaxy LOS distribution in redshift 
space is the measurement error of the SDSS galaxy red- 
shifts. We show in Appendix A an accurate modelling of the 
redshift errors from repeat observations of galaxy spectra 
in the SDSS. We find that the additional velocity contribu¬ 
tion introduced by the redshift errors is best modelled by 
a Gaussian-convolved Laplace distribution. We adopt two 
different redshift error models for the luminous and faint 
galaxies (see details in Appendix A). The typical la redshift 
error is about 10-15 kms“^. The redshift errors (following 
the Gaussian-convolved Laplace distribution) are built into 
our HOD model. 

Following G15, we apply a Markov Chain Monte Carlo 
method to explore the HOD parameter space. The likeli¬ 
hood oc exp(—x^/2) for a given set of HOD parameters is 
determined by the , contributed by the projected 2PCF 
WpRp), the redshift-space multipoles ^o{s), ^ 2 ( 3 ) and ^ 4 ( 3 ), 
and the observed galaxy number density Ug, 


x^ = (4-rrc-i(4-r) + (s) 


where C is the full error covariance matrix and the data vec¬ 
tor ^ = [lOp,^ 0 ,^ 2 ,^ 4 ]. The quantity with (without) a su¬ 
perscript is the one from the measurement (model). The 
covariance matrix is determined from 400 jackknife samples 
as mentioned above (Zehavi et al. 2011; Guo et al. 2013). 
We apply a mean correction for the bias effect in inverting 
the covariance matrix, as described in Hartlap et al. (2007) 
(see also Percival et al. 2014). We also apply a volume cor¬ 
rection of 1 -k V/Vsim to the covariance matrix to account 
for the model uncertainty caused by the the finite volume 
(Kim = 1 h~^GpC) of the MDPL simulation (G15). The er¬ 
ror arig on the number density is determined from the vari¬ 
ation of rig in the different jackknife samples. The volume V 
and mean number density rig of each luminosity threshold 
sample are listed in Table 1. 
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Figure 3. HOD fittings to the four sets of 2PCF measurements in Fig. 2. The measurements from the data are shown by open circles, 
while the HOD model fits are displayed as lines. For clarity, the measurements of different luminosity threshold samples are separated 
by 0.2dex (for Wp) or 20( /i“^Mpc)^ (for s^5o,2,4) per sample, starting from the Mr < —18 sample. In the top left panel of Wp{rp), the 
measurements at rp < 0.1 h~^Mpc (below the fibre-collision scale) are also shown, which are not included in the HOD fittings, while the 
dashed lines are the predictions from the best-fitting HOD models. 


4 RESULTS 

4.1 Fitting Results and the Mean Occupation 
Function 

Fig. 3 shows the measurements and best-fitting HOD models 
for the four sets of 2PCFs as in Fig. 2. For clarify, offsets 
are applied to both the data points and bestfitting curves, 
0.2dex for Wp and 20(/i“^Mpc)^ for s^^o, 2,4 for each galaxy 
sample. As is evident, our HOD model leads to remarkably 
good fits to all the luminosity threshold samples, for both 
the projected and redshift-space 2PCFs. We choose to fit 
2PCFs on scales above 0.1 h“^Mpc to reduce any possible 
systematic effect in the fibre-collision correction (which is 
small according to Guo et al. 2012). With our bestfitting 
HOD model, we can predict the 2PCFs on smaller scales. 
The dotted curves and the filled circles in the top left panel 
of Fig. 3 show the prediction and the measurements on scales 
below 0.1/i“^Mpc. The best-fitting models reproduce well 
the measurements for all luminosity threshold samples down 


to Tp = 0.02 h ^Mpc, suggesting both an accurate HOD 
model and a robust fibre-collision correction. 

The best-fitting HOD parameters are listed in Table 2 
for the different luminosity threshold samples. The x^/dof of 
the fittings confirm the adequacy of the model in fitting the 
data. The derived parameters Mi and /sat are also displayed, 
together with the characteristic central and satellite galaxy 
velocity dispersions CTc = acU'„(Mmin) and CTs = asU„(Mi) 
(see § 4.2 for more details). 

In our modelling, we choose to fit the projected 2PCF 
and the redshift-space 2PCF multipole ^ 0 , 2 , 4 , not the 3D 
redshift-space 2PCF 5(rp,ra-) directly. The reason is the 
large dimension of ^(rp,ra-) (e.g. 240 data points per sam¬ 
ple for 12 Tp bins and 20 r^- bins), which makes it difficult 
to estimate a robust covariance matrix. However, with the 
best-fitting models, we can predict C(^pU'n-) and compare 
to the measurements as a cross-check. Such a comparison is 
shown in Fig. 4 for three representative luminosity-threshold 
samples of Mr < —19, —20 and —21. Two main RSD ef- 
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Table 2. Best-fitting HOD parameters and derived parameters for the luminosity threshold samples 


M” 


-18.0 


-18.5 


-19.0 


-19.5 


- 20.0 


-20.5 


- 21.0 


-21.5 



Figure 4. Comparisons between the measured redshift-space 3D 2PCF ^(rp,r-^) (color scales) and the prediction from the best-fitting 
HOD model (black solid curves), for three representative luminosity-threshold samples. The data and models have the same contour 
levels of 0.3, 0.5, 1,2, 5, 10, and 20. 


fects show up in ^(rp,r.n.). On large scales, the galaxy in¬ 
fall towards overdense regions as well as the streaming of 
galaxies out of underdense regions compresses the contours 
along the LOS direction, known as the Kaiser squashing ef¬ 
fect (Kaiser 1987; Hamilton 1992). On small scales, the ran¬ 
dom motions of galaxies in virialized structures cause the 
contours to appear stretched along the LOS direc¬ 
tion, causing the FOG effect (Jackson 1972; Huchra 1988). 
The best-fitting HOD models reproduce the two features 
and the overall ^(rp, r,r) measurements very well. In partic¬ 
ular, the agreement between the measured and predicted 3D 
2PCFs on small scales is remarkably good. 

On large scales, there appears to be slight deviations of 
the predictions from the measurements, which are not signif¬ 
icant, given the measurement errors. To see this and to have 
a comparison for all samples, in Fig. 5 we compare the 3D 
2PCF contours for both the best-fitting predictions (red dot¬ 
ted curves) and the meausurements (blue solid curves). To 
illustrate the uncertainties on the measurements of CRpt 


we show in each panel with black dotted curves the contours 
from ±1(7 of the ^(rp,r.,r) measurements for the outmost 
level {^{rp,rTr) = 0.3). The model fits of the Mr < —18.5 
and Mr < —19 samples seem to be out of the la range of 
the measurements on large scales of Tp > 10/i~^Mpc, which 
is consistent with the over-prediction of the model on large 
scales, as seen in Wp and in Fig. 3. However, such devia¬ 
tions are not significant given the highly correlated covari¬ 
ance matrix elements on these scales. The largest difference 
in the FoG feature is seen in the Mr < —18.0 sample, which 
is in fact not significant given the large error bars in the mea¬ 
surement for this sample. Overall the model successfully re¬ 
produces the luminosity dependent ^(rp,r.,r) measurements. 
As will be discussed in § 4.2, velocity bias is needed for the 
model to fit the redshift-space clustering. 

The mean halo occupation functions for the different 
luminosity threshold galaxy samples are presented in Fig. 6. 
The most significant trend is that the host halo mass scale 
increases with the galaxy luminosity, as expected from HOD 
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Figure 5. Comparisons between the predicted 3D 2PCF ^(rp,r 7 r) from the best-fitting HOD models (red dotted curves) and the 
measurements (blue solid curves) for different luminosity-threshold samples. Contour levels are set to ^(rp, r-^) = [0.3, 0.5,1, 2, 5,10, 20]. 
The black dotted curves in each panel are contours at the level of $(rp,r 7 r) = 0.3 with $(rp,r 7 r) shifting by dilcr measurement error. The 
measurements of ^(rp,r 7 r) in each panel (as well as the model curves) are smoothed with a Gaussian kernel to reduce the noise on large 
scales. 



M (h-iMj 


of 

Figure 6. Mean halo occupation functions of the best-fitting 
models for different luminosity threshold samples. 


modelling of projected 2PCF Wp{rp) from the same SDSS 
data by Zehavi et al. (2011) (see also Zehavi et al. 2005; 
Zheng et al. 2007). Compared to the HOD model used in 
Zehavi et al. (2011), our model in this paper adopts different 
cosmological parameters and halo definition. Furthermore, 
we shift to a simulation-based model, rather than an analytic 
model. Finally, we jointly fit the projected 2PCF Wp and the 
redshift-space 2PCFs Co, 2 , 4 , while Zehavi et al. (2011) only 
fit Wp. Accounting for these differences, our results are in 
good agreement with those in Zehavi et al. (2011). We note 
that the nncertainties in many HOD parameters (and the 
derived satellite fraction) from the modelling in this paper 


appear to be larger than those in Zehavi et al. (2011) from 
modelling Wp only. This can be attributed to the differences 
in the models. The accuracy of the measured small-scale 
data points exceed that of the analytic HOD model used in 
Zehavi et al. (2011), which may be the reason of their large 
X^/dof (2-3 for some cases). As a consequence, the uncer¬ 
tainties in the HOD parameters can be artificially underes¬ 
timated. The simulation-based model used in this paper is a 
more accurate model, leading to good values of x^/dof and 
improved error estimates in the parameters. 

We also perform uip-only fit with the simulation-based 
model and compare to the results from fitting both Wp and 
Co, 2 , 4 - We find that redshift-space 2PCFs help tighten the 
constraints on the HOD parameters. As an example, we 
show in Fig. 7 the comparison of the constraints on Mmin and 
uiog M from fitting Wp only (blue contours) and jointly fitting 
Wp and Co, 2,4 (red contours). We set a prior of uiogM < 1-5 
when fitting the data to have a reasonable value of the scat¬ 
ter. Clearly, a substantial improvement with the redshift- 
space 2PCFs is to narrow down the range of awgM, espe¬ 
cially for less luminous samples (with the Mr < —18 as an 
exception, which has a tighter Mmin). 

Even though redshift-space 2PCFs help tighten the con¬ 
straints on (JiogM, we note that for faint galaxy samples the 
cutoff profile in the mean central occupation function is still 
not well constrained, as indicated by the large errors (Ta¬ 
ble 2 and Fig. 7). It is consistent with a sharp cutoff at 
Mmin, and in Fig. 6 we choose to plot the best-fitting mod¬ 
els with CTiogM — 0 for these samples. The constraints on 
Mmin and (JiogM mainly come from the galaxy bias (large 
scale 2PCF amplitude) and the galaxy number density. The 
galaxy bias is mainly determined by haloes around Mmin- 
For faint samples, Mmin is in the range that halo bias is in¬ 
sensitive to halo mass. As a consequence, the galaxy bias is 
insensitive to the way of populating galaxies into haloes of 
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Figure 7. Comparisons between the constraints on the HOD parameters Mmin a-iid criog M from fitting Wp (blue contours) and jointly 
fitting Wp and ^ 0 , 2,4 (red contours). 


different masses around Mmin, i.e. insensitive to the change 
in (JiogM- A change in criogM can be easily compensated 
by a slight change in Mmin to maintain the galaxy num¬ 
ber density. Therefore, the cutoff profiles for faint samples 
are not well constrained. Conversely, ciogM is much better 
constrained for the luminous samples as a result of the steep 
dependence of halo bias and halo mass function on halo mass 
toward the high-mass end. 


Fig. 8 shows the dependence of the characteristic mass 
scales (Mmin for central galaxies and Mi for satellite galax¬ 
ies) and the satellite fraction /sat on the sample number 
density Ug. The dependence of any of those parameters 
on rig roughly follows a power-law form. As pointed out 
in Guo et al. (2014), the dependence of Mmin on the num¬ 
ber density largely comes from the nearly power-law form 
of the halo mass function over a large mass range. The 
mass Mmin is mostly determined by matching the halo num¬ 
ber density with the galaxy number density, modulated by 
criogM. There is a trend that the ratio Mi/Mmin decreases 
as the sample number density decreases (or the sample lu¬ 
minosity increases), consistent with those found in litera¬ 
ture (e.g. Zehavi et al. 2005; Skibba et al. 2007; Zheng et al. 
2009; Zehavi et al. 2011; Guo et al. 2014; McCracken et al. 
2015; Skibba et al. 2015). This is a manifestation of the halo 
mass dependent competition between accretion of galaxies 
into haloes and destruction of galaxies inside haloes. From 
the bottom panel of Fig. 8, we see that fainter galaxies are 
more likely to be satellite galaxies in massive haloes. The 
satellite fraction follows /sat — 0.1[hg/(10“®/i®Mpc“®)]^^^, 
as also shown in other surveys (Guo et al. 2014), which can 
be used to estimate the satellite fraction given the number 
density of a threshold galaxy sample. Note that the values 
of satellite fraction are slightly lower than those inferred in 
Zehavi et al. (2011), which can be mostly attributed to the 
difference in the definitions (hence the sizes) of haloes. 



0.1 10-2 lQ-3 10 -'' 

rig (h^ Mpc-2) 

Figure 8. Top: characteristic mass scales (Mmin and Mi) of 
haloes hosting central galaxies and satellites as a function of the 
sample number density Ug. Bottom: satellite fraction as a func¬ 
tion of the sample number density. The shaded area around each 
curve shows the Icr uncertainties in the parameter. 
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4.2 Galaxy Velocity Bias 

The constraints on the galaxy velocity bias parameters are 
shown in Fig. 9, including the 68 per cent and 98 per cent 
contours in the oic-Us plane and the marginalized distribu¬ 
tions of Oc and as for each sample. In terms of the tightness 
in the central galaxy velocity bias Oc constraints, the lumi¬ 
nous and faint samples show a dichotomy. 

For the luminous samples (more luminous than Mr = 

— 19.5), both the central and satellite velocity bias param¬ 
eters are well constrained, as shown in the top left panel. 
The case without any galaxy velocity bias (i.e. Qc = 0 and 
Os = 1) is far beyond the 95 per cent contours of all lu¬ 
minous samples. That is, galaxy velocity bias is required to 
reproduce the redshift-space clustering in the local universe 
for luminous samples. 

The central velocity bias parameter a^ for luminous 
samples is about 0.3 (top-left and bottom-left panels in 
Fig. 9). It shows no significant dependence on galaxy lu¬ 
minosity. The existence of the central velocity bias implies 
that these luminous central galaxies are not at rest at the 
halo centres with respect to the bulk motion of the haloes. 
This reflects the mutual (non)relaxation status of central 
galaxies and host haloes, which are related to the merger 
history of the galaxies and the formation history of haloes 
(see more detailed discussions in G15). The value of Oc in¬ 
ferred from our modelling of the redshift-space clustering is 
in agreement with the estimates from galaxy group cata¬ 
logues in the SDSS (e.g. van den Bosch et al. 2005), which 
uses the mean velocity of satellites as a proxy for the halo 
velocity. 

The constraints on the central velocity bias parameter 
Qc are loose for the three faint samples (with threshold lumi¬ 
nosity fainter than Mr = —19.0), as seen from the contours 
in the top-right panel and the corresponding curves in the 
bottom-left panel of Fig. 9. For the Mr < —19 sample, ac is 
different from zero only at the 2.5cr level (see Table 2). For 
the other two fainter samples, a^ is consistent with zero. The 
loose constraints can be partly attributed to the relatively 
large uncertainty in the clustering measurements and in the 
jackknife covariance matrix estimate for the faint samples, 
as the sample volumes are substantially smaller than those 
of the luminous samples (see Fig. 1 and Table 1), especially 
on large scales where a^ is mostly constrained (see Fig.6 of 
G15). The three faint samples have volumes that are 10%, 
24%, and 43% that of the Mr < —19.5 sample, which has 
the smallest volume among the luminous samples. The other 
possible cause of the loose constraints can be the redshift 
errors. As the central velocity bias, in terms of velocity dis¬ 
persion (see below), approaches or drops below the level of 
redshift errors (about 13 kms“^; see Appendix A), the sen¬ 
sitivity of RSD to the central velocity bias is reduced, likely 
the case for the faint samples. 

For the satellite velocity bias, all samples show good 
constraints. Sample volume becomes less important here, 
since the constraints mainly come from the small-scale FoG 
effect (see Fig.6 of G15), where the uncertainties in the mea¬ 
surements are small. However, the faintest sample {Mr < 

— 18) may still be affected by the small sample size and 
the noisy covariance matrix estimate, and we should inter¬ 
pret the satellite velocity bias with caution. If we neglect 
the Mr < —18 sample, the satellite velocity bias as con¬ 


straints show a more or less continuous trend with lumi¬ 
nosity, i.e. higher values of as for more luminous sample. 
Probably more appropriately, the satellite velocity bias con¬ 
straints can be divided into two groups. For the two most 
luminous samples (more luminous than L*, corresponding to 
Mr = —20.44; Blanton et al. 2003b), as is consistent with 
unity. That is, the motion of the luminous satellites closely 
follows that of the dark matter. For samples with threshold 
luminosity fainter than L*, as is about 0.8 “ 0.85. That is, 
for those samples, satellites moves more slowly than dark 
matter particles. 

In a steady state, the spatial distribution and velocity 
distribution of satellite galaxies inside haloes are related to 
each other. In our modelling, we draw random dark matter 
particles for the position of satellites. That is, we implic¬ 
itly assumed that the spatial distribution of satellites fol¬ 
lows that of the dark matter, which is well described by the 
Navarro-Frenk-White (NFW) profile (Navarro et al. 1997). 
For the two luminous samples, as is around unity, i.e., their 
satellite velocity distribution is consistent with that of the 
dark matter. Therefore, based on the constraints we infer, an 
NFW profile and associated velocity distribution for satel¬ 
lites are able to explain the redshift-space clustering for the 
luminous sample. More modelling efforts are needed to see 
whether other profiles and the corresponding velocity dis¬ 
tributions are preferred or not (see the tests in G15). For 
the other, faint samples, the inferred as (~0.8 - 0.85) dif¬ 
fers substantially from unity, inconsistent with the value for 
the NFW profile. The result alone suggests that the spatial 
distribution of faint satellites should deviate from the NFW 
profile. Our current model, however, is not able to provide 
more information on how significant a deviation it needs to 
be. For an improved model, one can consider to parameterize 
the spatial profile of satellites and solve for the correspond¬ 
ing velocity distribution in a self-consistent manner, which 
is beyond the scope of this paper. 

Watson et al. (2010) and Watson et al. (2012) analyzed 
the small-scale (down to ~0.01 /i“^Mpc) clustering {w-p) of 
the SDSS Main galaxy sample and luminous red galaxies. 
They found that the spatial distribution of faint satellite 
galaxies (below Mr = —20) is consistent with the NFW 
profile, while that of bright satellite galaxies deviates from 
the NFW profile (with a steeper Inner profile). These seem 
opposite to what we find from modelling the redshift-space 
clustering. An Improved model Is necessary to constrain the 
range of profiles allowed by the redshift-space clustering data 
and to see whether this apparent difference Is significant. 

In the left panel of Fig. 10, we summarize the con¬ 
straints on the velocity bias parameters a^ and Qs as a func¬ 
tion of galaxy number density (more luminous galaxy sam¬ 
ples have lower number densities). For galaxies Mr < —19.5 
and brighter, there Is no significant dependence of a^ on 
number density. The three faint galaxy samples (with the 
highest number densities) seem to have lower values of Oc, 
consistent with zero, but with large error bars. Since the host 
halo mass Increases as the galaxy luminosity Increases (or as 
the number density decreases), the trend also Indicates the 
weak dependence of a^ on the host halo mass. With the 
faintest sample excluded (smallest sample volume), the de¬ 
pendence of satellite velocity bias as on the luminosity for 
faint galaxies Is also weak, while as shows a clear Increase 
with luminosity for Mr < — 21. 
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Figure 9. Top panels: marginalized probability distributions of central and satellite galaxy velocity bias parameters for the luminous (left) 
and faint (right) galaxy samples. The contours show the 68 and 95 per cent confidence levels for the two parameters. Bottom panels: ID 
probability distribution of the central (left) and satellite (right) galaxy velocity bias parameters for various luminosity-threshold samples. 


To interpret the velocity bias results, the more mean¬ 
ingful physical quantities are the velocity dispersions of 
central and satellite galaxies inside haloes, denoted as 
cTc and (Ts, respectively. Given the velocity bias param¬ 
eters, velocity dispersions depend on halo mass, and we 
choose to evaluate typical values in representative haloes. 
For central galaxies, the velocity bias constraint is mainly 
contributed from haloes around Mmin, and we compute 
the typical central galaxy velocity dispersion as Uc = 
ac(T„(Minin). For satellite galaxies, the effect of the ve¬ 
locity bias on clustering comes from haloes around Mi, 
and the typical satellite velocity dispersion is computed 
as (Ts = asCJv{M\). The right panel of Fig. 10 shows the 
dependencies of these typical velocity dispersions on sam¬ 
ple number density, which roughly follow power-law rela¬ 
tions, Oc — 55kms“^[hg/(10“®h^Mpc“®)]“°'^® and (Ts — 
437kms"^[h3/(10"^h®Mpc"^)]"^/®. 

The existence of galaxy velocity bias reflects the dynam¬ 
ical evolution of galaxies inside haloes. For example, infalling 
satellites experience tidal striping and dynamical friction in 
the dark matter haloes, affecting its velocity distribution. 


For central galaxies, the existence of velocity bias indicates 
that central galaxy and the host haloes are not mutually 
relaxed. A likely cause can be the halo mergers and the sub¬ 
sequent galaxy mergers. Our results can be used to assess 
the dependence of the degree of relaxation after mergers on 
halo mass (or galaxy number density). For mergers of haloes 
of similar mass, the mean pairwise infall velocity ui 2 on large 
scales is proportional to the bias factor (Sheth et al. 2001; 
Zhang & Jing 2004). With the luminosity-dependent bias 
factor in Zehavi et al. (2011), we find that the bias factor ap¬ 
proximately scales as , which means that the pairwise 

infall velocity vi 2 oc before merger. The central ve¬ 

locity dispersion constrained from the RSD iotc oc n“°'^®) is 
much steeper than the infall velocity. We therefore conclude 
that central galaxies in lower mass haloes are more relaxed 
with respect to the host haloes, compared with their coun¬ 
terparts in more massive haloes, consistent with an overall 
earlier formation and thus more time for relaxation of the 
lower mass haloes. 
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Figure 10. Velocity bias and typical velocity dispersion. Left: central and satellite galaxy velocity bias as a function of the sample 
number density rig. Right: dependence of the typical central and satellite galaxy velocity dispersions, ctc and as, on the sample number 
density. The shaded area show the 68 per cent and 95 per cent confidence levels of the parameters. The solid lines with circles are the 
best-fitting models. 



Figure 11. Dependence of central and satellite galaxy velocity 
bias on the cosmological parameter Qm- The central and satel¬ 
lite velocity bias are displayed in circles of different colours. The 
shaded area show the 68 and 95 per cent confidence levels of the 
parameters. 

4.3 Dependence of Velocity Bias on Cosmology 

The velocity bias constraints we infer rely on the Mul- 
tiDark simulation with the assumed cosmology. The cos¬ 
mological parameters are close to the results from Planck 
(Planck Collaboration et al. 2014). It is useful to see 
whether the velocity bias constraints are robust against rea¬ 
sonable change in cosmology. Here, we limit our investiga¬ 
tion to the change in Instead of using simulations with 
different parameters, we use the appropriate scaling re¬ 
lations with the MDPL simulation for the corresponding 
change in the dark matter halo properties when varying 

According to Zheng et al. (2002) and Tinker et al. 


(2006), if two simulations have the identical initial matter 
fluctuation spectrum (including amplitude and phase) but 
different values of ilm, there exists a simple relation between 
simulation outputs at a given linear growth factor G. There 
are correspondences between haloes in the two simulations, 
and the corresponding haloes have the same radius with the 
mass scaling with 12^. Halo velocity scales with the growth 
rate / = dlnG/dlna, with a the scale factor. The internal 
velocity dispersion inside haloes scales with We there¬ 
fore modify our default simulation output by scaling the halo 
mass M, halo velocity Vh, dark matter particle velocity Vp, 
and the ID velocity dispersion (j„ of dark matter particles 
in haloes in the following way. 


M 

= Mo, 

(9) 

Vh 

= {f/fo)Vh,0, 

(10) 

- Vh 


(11) 

(J V 

— ^ O'v,0f 

(12) 


where the symbols with subscript ‘0’ denote the values in 
the fiducial MDPL simulation. The halo mass function also 
changes accordingly in each scaled simulation. We build 
eleven scaled simulation catalogues, varying Dm from 0.24 
to 0.34 with a step of Aflm = 0.01. We then apply our HOD 
model to the Mr < —21 sample based on the eight scaled 
simulations to constrain the velocity bias parameters. 

We show in Fig. 11 the dependence of the velocity bias 
parameter constraints on Dm. Both the central and satellite 
velocity biases decrease with increasing Dm, as expected. 
Increasing Dm leads to higher halo-halo and internal halo 
velocity dispersions, and the velocity bias parameters de¬ 
crease to compensate such a change to match the redshift- 
space clustering. For the range of Dm considered here, both 
velocity bias parameters are well constrained. In particular, 
the central velocity bias differs from zero for all the cases, 
which shows that our constraints are robust against reason¬ 
able changes in cosmology. 

We note that the luminosity-threshold samples used 
in this paper are constructed over different redshift ranges 
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according to the range for which they are volume limited 
(see Fig. 1 and Table 1). If we were to apply the cosmol¬ 
ogy change and compare the results among different sam¬ 
ples, the difference in the sample mean redshifts needs to 
be accounted for, which leads to small effective cosmology 
changes. As mentioned in § 3, our test with the model built 
on the 2 = 0.1 simulation output shows that the inferred 
HOD parameters are within the uncertainties of those from 
the default model built on the 2 = 0 output. So the effect 
of cosmology change from the difference in mean redshift 
is small. Furthermore, here we only focus on one sample 
and aim to see the sensitivity of velocity bias constraints to 
cosmology, and the effective cosmology change among dif¬ 
ferent samples becomes irrelevant. With one sample, there 
is another effect related to the cosmology change. In prin¬ 
ciple, when varying the cosmological parameter Q.m, we 
need to rescale or remeasure the 3D 2PCFs and then per¬ 
form the modelling. Otherwise, the Alcock-Paczynski effect 
(Alcock & Paczynski 1979) would introduce an additional 
distortion in the 2PCFs. However, since the galaxy samples 
are local, with 2 ^ 0 , the comoving distance is insensitive to 

and the effect is tiny. For example, at the mean redshift 
of the Mr < —21 sample, 2 ~ 0.12, the comoving distance 
changes by less than 1% in the whole range of Qm we study 
here, which has little effect on the 2PCF measurements. We 
therefore do not rescale or remeasure the 2PCFs changing 
the cosmology. 

Relevant to the RSD effect, a more interesting change 
in cosmology is in the combination of the growth rate and 
fluctuation amplitude, fas. Because of the high precision 
measurements, the small-scale redshift-space clustering can 
help tighten the constraints on fag (e.g. Reid et al. 2014), 
providing potentially stringent tests to the ACDM cosmol¬ 
ogy and theory of gravity. We reserve such an investigation 
on the fag constraints with the SDSS Main galaxies for a 
future work. 


5 SUMMARY AND DISCUSSIONS 

We measure the projected and redshift-space 2PCFs for 
volume-limited, luminosity-threshold samples of SDSS Main 
galaxies, on small to intermediate scales (0.1-25/i“^Mpc). 
The measurements are interpreted within the HOD frame¬ 
work to infer the relation between galaxies and dark matter 
haloes. In particular, the RSD effects in the redshift-space 
2PCFs enable us to constrain the kinematics of central and 
satellite galaxies inside dark matter haloes and infer the dif¬ 
ference between the motions of galaxies and dark matter. 

It is the first time that the redshift-space clustering 
of local galaxies is accurately measured on scales as small 
as ~0.1/i~^Mpc. The measurements become possible with 
the accurate fibre-collision correction method developed in 
Guo, Zehavi & Zheng (2012), which makes use of the re¬ 
solved collided galaxy pairs in tile overlap regions to re¬ 
cover the small scale clustering. Previous measurements (e.g. 
Hawkins et al. 2003; Zehavi et al. 2005) rely on either an¬ 
gular or nearest neighbour fibre-collision corrections, which 
results in systematics at the level of the data precision 
(Guo et al. 2012). With our measurements, we find that 
both the projected and redshift-space 2PCFs show a clear 
dependence on galaxy luminosity, generally with a higher 


clustering amplitude for more luminous galaxies on both 
small and large scales. The dependence on luminosity be¬ 
comes stronger for galaxies above ~L*. The overall trend 
is consistent with previous results based on the projected 
2PGFs (e.g. Zehavi et al. 2011). 

To interpret the measurements, similar to G15, we re¬ 
sort to an accurate HOD model based on a high-resolution 
A-body simulation. In addition to the mean halo occupa¬ 
tion function, the model also parameterizes the central and 
satellite galaxy velocity bias. For the first time, a halo-based 
model is applied to model the measured luminosity depen¬ 
dent small- and intermediate-scale redshift-space clustering 
of local galaxies. Previous studies usually focus on relatively 
large scales and adopt a streaming model (e.g. Peacock et al. 
2001; Hawkins et al. 2003; Bel et al. 2014; Hewlett et al. 
2015). The commonly inferred quantities include the linear 
redshift distortion parameter and the mean pairwise velocity 
dispersion of galaxies (e.g Cabre & Gaztafiaga 2009) or its 
scale dependence (e.g. Li et al. 2007). Our halo-based model, 
as applied in G15 (see Reid et al. 2014 for a similar model), 
makes use of the kinematic information of haloes and param¬ 
eterizes the galaxy velocity distribution on top of it, rather 
than an overall mean velocity dispersion. It allows us to con¬ 
strain the occupation and kinematic distribution of galaxies 
at the level of dark matter haloes, a more informative extrac¬ 
tion from galaxy redshift-space clustering data. We hnd that 
the model is able to successfully reproduce the observed pro¬ 
jected 2PCF Wp{rp), the redshift-space multipole moments 
^o(s), ^ 2 (s) and ^ 4 (s), and the 3D 2PCF f{rp,rr), on all 
scales for all SDSS luminosity-threshold samples (Fig. 3 and 
Fig. 5). 

Consistent with previous work that only model the pro¬ 
jected 2PGF (e.g. Zheng et al. 2007; Zehavi et al. 2011), we 
find that the clustering trend with luminosity can be ex¬ 
plained by the fact that more luminous galaxies reside in 
more massive haloes. The characteristic halo masses, Mmin 
(for central galaxies) and Mi (for satellite galaxies), increase 
with increasing luminosity threshold. The satellite fraction, 
/sat, drops as the luminosity threshold increases. The depen¬ 
dence of Mmin, Ml, and /sat on the sample number density 
(which is directly related to the sample luminosity thresh¬ 
old) can be well described by power-law relations. Gompared 
to the Wp-only modelling results, the redshift-space 2PGFs 
help tighten the constraints on the HOD parameters. How¬ 
ever, we find that for the faint galaxy samples (with lumi¬ 
nosity threshold below L*), the cutoff prohle (characterized 
by the parameter uiogM) in the mean occupation function 
of central galaxies is still loosely constrained. 

Besides the above results, the brand-new outcomes from 
our modelling are the constraints on galaxy kinematics in¬ 
side haloes, coming from the RSD effects on both small 
and large scales. The redshift-space clustering data require 
the existence of a non-zero central galaxy velocity bias of 
about 0.3 for luminous samples (with threshold luminosity 
Mr < —19.5), while for faint samples the central velocity 
bias parameters are loosely constrained but consistent with 
the above value. That is, in the rest-frame (centre-of-mass 
frame) of a halo, the central galaxy on average moves at a 
speed about 30% that of dark matter particles. The cen¬ 
tral galaxy velocity bias in our model is in agreement with 
the estimates from the galaxy group catalogues in the SDSS 
(e.g. van den Bosch et al. 2005). This mutual non-relaxation 
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between central galaxies and dark matter haloes can re¬ 
sult from mergers and dynamical evolution of galaxies and 
haloes. Converting the value to physical speed and compar¬ 
ing to the typical infall velocity before merging, our results 
imply that galaxies in lower mass haloes are more relaxed 
with respect to the host haloes, consistent with their earlier 
formation time. Further theoretical investigation of such an 
evolution paradigm using high-resolution cosmological hy¬ 
drodynamic galaxy formation simulations will be pursued 
in future work. 

For satellite galaxies, we find that the two most lu¬ 
minous samples have satellite velocity bias consistent with 
unity, which means that they follow closely the motion of 
dark matter. The satellite velocity bias ois ~ 0.85 for fainter 
samples implies that fainter satellites move more slowly than 
dark matter. If satellite motion is in a steady state, the re¬ 
sult suggests that the spatial distribution profile of faint 
satellites should differ from that of dark matter. More in¬ 
formative constraints need an improved model that allow 
for a self-consistent treatment of the spatial and kinematic 
distributions of satellites, which would also lead to tighter 
constraints than using only projected spatial clustering mea¬ 
surements (e.g. Watson et al. 2012; Wang et al. 2014b). An 
improved model can also include the effect of halo assem¬ 
bly bias to see its influence on the HOD (e.g. Zentner et al. 
2014; but see Lin et al. 2015) and velocity bias. The exis¬ 
tence of satellite galaxy velocity bias affects any dynamical 
inference based on satellite velocity dispersions. For exam¬ 
ple, the effect needs to be taken into account when using 
the velocity dispersion of the galaxy cluster members to es¬ 
timate the halo mass of galaxy clusters (Goto et al. 2003; 
Old et al. 2013). 

In G15, velocity bias is inferred for a sample of z ^ 
0.5 luminous galaxies {rig = 2.19 x 10“'^h®Mpc“®). Con¬ 
verted to the same velocity bias definition used in this pa¬ 
per (i.e. in the centre-of-mass frame of haloes), their re¬ 
sult on the velocity dispersion for central galaxies in haloes 
of ~ is about 96 ± 8kms“^. On average, 

these haloes evolve to ~ ]^q13-55 /j-IjvJq haloes at z~ 0 (e.g. 
Zhao et al. 2009), around Mmin of the Mr < —21.5 sample 
{ug = 2.86 X 10~^/i^Mpc“^). We thus can have an approx¬ 
imate connection between the two samples at a ~ 0.5 and 
z ^ 0 from the halo evolution, following the same spirit as in 
Zheng et al. (2007), which is also supported by the number 
density comparison. The central galaxy velocity dispersion 
for the Mr < —21.5 is about 99± 15 kms“^. Such a prelimi¬ 
nary analysis shows that the velocity dispersion for luminous 
central galaxies has not evolved much since 2 ~ 0.5 (in fact, 
it has marginally increased, but only at a ~ 0.3cr level). We 
can assume a simple circular motion of the central galaxy in 
the inner NFW halo to study the implication of the result. 
With the 2 ~ 0.5 velocity bias result, the radius of the orbit 
can be inferred to be about 0.4% of the virial radius of the 
halo (G15). The corresponding dynamical friction time scale 
is estimated to be ~0.1 Myr, much shorter than the 3.7 Gyr 
time span from 2 ~ 0.5 to 2 ~ 0.1. Therefore, a substantial 
central velocity bias at 2 ~ 0 indicates that these luminous 
central galaxies and their host haloes (or the cores of haloes 
and the rest of the haloes) may have been constantly dis¬ 
turbed by galaxy and halo mergers. We expect that velocity 
bias of galaxies inferred from different redshifts can help 


study the dynamical evolution of galaxies and haloes and 
test galaxy formation theory. 

The galaxy peculiar velocity field is directly related 
to the growth of structures in the universe. RSD effects 
can be used to constrain cosmology, especially the growth 
rate (i.e. the fas parameter) to test the theory of gravity 
(e.g. Guzzo et al. 2008; Percival & White 2009). The statis¬ 
tical power of small-scale RSD measurements have the po¬ 
tential to greatly tighten the constraints (e.g. Tinker et al. 
2006; Tinker 2007; Reid et al. 2014). The effect of velocity 
bias needs to be taken into account when using small-scale 
redshift-space clustering to constrain cosmology. We leave 
an investigation on cosmological applications based on our 
measurements and model for a future work. 
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APPENDIX A: THE REDSHIFT ERROR 
DISTRIBUTION 

In order to develop an accurate HOD model to apply to the 
redshift-space clustering measurements in SDSS, we need 
to carefully account for the distribution of redshift measure¬ 
ment errors of galaxies. The effect of redshift errors is to add 
apparent peculiar velocity dispersion to galaxies, which, if 
not accounted for, introduces apparent velocity bias compo¬ 
nent. This is especially important for faint galaxies, because 
the potentially small central and satellite velocities inside 
haloes make the constraints more vulnerable to redshift er¬ 
rors. 

To properly investigate the redshift errors in the SDSS 
Main galaxies, we make use of all the galaxies with repeat 
spectra for each luminosity-threshold sample to derive the 
distribution of redshift measurement errors. In particular. 
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Figure Al. Redshift error distribution of SDSS Main galaxies with different luminosity thresholds. In each panel, the histogram 
is obtained from galaxies with repeat spectroscopic observations, and the red curve is the best-fitting Gaussian-convolved Laplace 
distribution. See the text. 


for a galaxy with n observations and hence n redshift mea¬ 
surements, we construct the estimator of the velocity error 
in each measurement as 


Avi = 


'c{zi - z) 


n — 1 


(Al) 


where Zi is the i-th redshift measurement for this galaxy 
{i = 1,..., n) and z is the mean of the n measurements, and 
the factor sjnjiji — 1) makes the estimator unbiased. 

The histograms in Fig. Al show the probability distri¬ 
bution of the redshift errors in terms of the peculiar velocity 
errors Av in four different luminosity-threshold samples. We 
find that the distribution has more extended tails than a 
naively assumed Gaussian distribution. The extended part 
follows closely the double exponential distribution or the 
Laplace distribution (appearing as straight lines in the plot), 
while the central part can be described by a Gaussian core 
(less sharply peaked than the Laplace distribution). In fact, 
the distribution can be remarkably fitted by a Gaussian- 
convolved Laplace distribution, shown as the red curves of 
Fig. Al. A random deviate Av for such a distribution can 
be obtained by the sum of two independent Gaussian and 
Laplace random numbers, i.e. Av = Auga,u +Auexp, with the 
distribution functions of Augau and Auexp as 


/gau(AUgau) 


/exp ( AUexp) 


■ exp 


■ exp 


-\/^o 


\f2o 


AUgau \ 
2cr|au / ’ 
v/2|AUexp| 

C^exp 


(A2) 

(A3) 


The parameters Ugau and Cexp are the standard deviations 
for the Gaussian and Laplace distributions, respectively. 

Why does the redshift error follow more closely to a 
Laplace distribution than a Gaussian distribution? Laplace 
distribution can be thought as the distribution of Gaussian 
random variables with mean zero and stochastic variance 
that has an exponential distribution. We speculate that for 
a given galaxy, the redshift error follows a Gaussian dis¬ 
tribution, but the variance of the distribution varies from 
galaxy to galaxy. We expect that for higher variance in red¬ 
shift errors, there are fewer galaxies, and that there is a 
minimum variance (e.g. from the spectral resolution). With 
such an expectation, the probability distribution of the vari¬ 
ance can be approximated as an exponential distribution 
with a cutoff towards lower value. The overall distribution 
of the redshift errors then follows a Laplace distribution 
with the central part modified by the lower cutoff on the 
variance. In such a scenario, the random variable Av can 
be obtained through the product of two independent ran¬ 
dom numbers, Av = au, where u follows the unit normal 
distribution and follows a modified exponential distri¬ 
bution of scale parameter with non-zero values above 
a threshold The two parameters CTe and at character¬ 
ize the tail and central part of the distribution. The above 
speculation is supported by the data. In Fig. A2, we show 
the histogram of the sample variance of the redshift er¬ 
ror distribution estimated from galaxies with n — 2 repeat 
spectra, with = 5D(Li[c(-2i — z)/{l -|- z)]^/{n — 1) for each 
galaxy. The dotted curve shows an exponential distribution 
of the variance of the Gaussian error distribution with 
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Figure A2. Distribution of sample variance of redshift errors. 
For each galaxy with n = 2 repeat spectra, the sample variance 
is computed as distribu¬ 

tion is consistent with the expected distribution (solid curve) 
for the following case: the redshift error of each galaxies follows a 
Gaussian distribution, while the variance of the Gaussian dis¬ 
tribution varies from galaxy to galaxy and this stochastic variance 
<7^ follows an exponential distribution with a cutoff (illustrated 
by the dotted curve, characterized by a scale parameter (t^ and 
a cutoff threshold ctj). The description provides an explanation 
for the redshift error distribution of SDSS galaxies (Fig. Al). See 
text for more details. 

(Tg = 200(kms“^)^ and cr^ = 20(kms~^)^. The expected dis¬ 
tribution of sample variance is shown as the solid curve. It 
is derived from convolving the dotted curve with the sample 
variance distribution at given [noting that (n — l)s^/cr^ 
follows a distribution with n — 1 degrees of freedom]. 
Clearly, the distribution from galaxies with repeat spec¬ 
tra is consistent with being from a sample of Gaussian errors 
with stochastic variance that follows a truncated exponen¬ 
tial distribution. 

In this paper, we adopt the Gaussian-convolved Laplace 
distribution to model the redshift error distribution. From 
fitting the histograms of redshift errors with such a distri¬ 
bution for different luminosity-threshold galaxy samples, we 
derive the parameters Ugau and Cexp. We find that it is suf¬ 
ficient to adopt two groups of parameters, for bright and 
faint galaxies, respectively. We use Ugau = 2kms“^ and 
Cexp = 12.5 kms“^ for the luminosity threshold samples of 
Mr < —18, —18.5 and —19, while the more luminous sam¬ 
ples have (Tgau = 5kms“^ and Cexp = 14.5 kms“^. We in¬ 
corporate the redshift errors into our model by adding shifts 
in the redshift-space positions of galaxies, following the cor¬ 
responding Gaussian-convolved Laplace distribution. 

This paper has been typeset from a ILjX/DTLjX file prepared by 
the author. 
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